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^ ■ Abstract 
00 

, Neutrino is a good probe of extra dimensions. Large mixing and the apparent 

lack of very complicated oscillation patterns may be an indication of large 



o 



couplings between the brane and a single bulk neutrino. A simple and re- 
alistic five-dimensional model of this kind is discussed. It requires a sterile 
in addition to three active neutrinos on the brane, all coupled strongly to 



Qh! one common bulk neutrino, but not directly among themselves. Mindful that 



sterile neutrinos are disfavored in the atmospheric and solar data, we demand 



induced mixing to occur among the active neutrinos, but not between the 
active and the sterile. The size R of the extra dimension is arbitrary in this 
model, otherwise it contains six parameters which can be used to fit the three 
neutrino masses and the three mixing angles. However, in the model those six 
parameters must be suitably ordered, so a successful fit is not guaranteed. It 
turns out that not only the data can be fitted, but as a result of the ordering, 
a natural connection between the smallness of the reactor angle #13 and the 
smallness of the mass-gap ratio AM^ olar / AM% tmospheric can be derived. 
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I. INTRODUCTION 



Neutrinos are different from the other fermions. They have small masses and large 
mixing. Small masses can be explained by the seasaw mechanism, or by the presence of 
large extra dimensions [ffl. Large mixing is usually obtained with a right texture, which in 
the presence of extra dimensions can be arranged through appropriate couplings between 
the three or four brane neutrinos and the three or four bulk neutrinos. In this paper, we 
discuss a new and economical mechanism to generate large mixing H, through the strong 
coupling of a single bulk neutrino to the brane neutrinos. 

We take the point of view that the main difference between neutrinos and other Standard 
Model (SM) fermions is due to the presence of extra dimensions. Being SM singlets, sterile 
neutrinos can roam in the bulk whereas other fermions cannot. If SM quantum numbers 
are trapped on the brane, perhaps flavor and generation numbers are trapped there as well, 
in which case there is no reason to require more than one neutrino in the bulk. With 
only one bulk neutrino to couple to, textures can no longer be arranged in the usual way, 
so the question is whether large and correct mixing between the brane neutrinos can be 
automatically induced by their couplings with this one common bulk neutrino. The answer 
is clearly no if the brane-bulk couplings are weak, but there is a chance when they are strong. 

To be definite we consider a simple model with one extra flat dimension of radius R, which 
we use as the inverse-energy unit in the problem. This model contains / brane neutrinos, 
each with its own Majorana mass rrii/R, and a single sterile bulk neutrino, coupled via 
Dirac masses di/R to the brane neutrinos. The brane neutrinos do not directly couple 
among themselves, and the bulk neutrino does not have a 5D mass. All phases are ignored 
in this simple model, and for technical simplicity we shall take all mi to be positive, and not 
equal to an integer nor a half integer. The / = 3 case of this model with equal couplings di 
was first proposed in Ref. 0. 

We will not ask in this paper the dynamical origin of the Majorana mass nor the coupling; 
we will simply take the 2f dimensionless parameters and di to be real and adjustable. 
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Unless that works phenomenologically, there is no point to speculate on the dynamical origin 
of the model. 

Models of this kind have been studied by many authors @-§|, but often without the 
Majorana masses, with / = 1, and/or in the weak coupling limit. When a single brane 
neutrino couples weakly to the bulk, oscillations occur between the brane neutrino and only 
a few low-lying Kaluza-Klein (KK) modes of the sterile bulk. Nevertheless, the participation 
of the few excited KK modes is important as they improve the poor energy dependence for 
a solar v e to oscillate into a sterile neutrino ||. To explain also atmospheric oscillation, it 
is necessary to introduce additional bulk neutrinos and arrange the texture of the couplings 
appropriately. 

When couplings are increased, more KK modes are involved, which causes a wild and 
complicated oscillation pattern to emerge ||. This could have been a smoking gun for 
the presence of an extra dimension, but unfortunately there is no hint in the experimental 
data to support it. The probability to convert an active neutrino into a sterile KK mode 
is also enhanced, a fact which might cause some discomfort because oscillation into sterile 
neutrinos is disfavored by experiments^ In the strong coupling limit when / = 1, the 
worrisome wild oscillations disappear, because the widths of the states are now larger than 
their separations, making the spectrum effectively a continuum. Leakage into the KK modes 
is now total, with a leakage rate increasing with the coupling |J, so very quickly the single 
brane neutrino simply disappear in sight. If / > 1, leakage is still present but is now divided 
among the / brane neutrinos. In addition to the continuum in the spectrum, there are now 
/ — 1 isolated eigenstates through which brane neutrinos with induced mixing can oscillate. 
For / = 3, there are 2 isolated eigenvalues and hence only one mass difference, making it 
impossible to explain solar and atmospheric oscillations simultaneously 0. 

In this paper we show that the strong coupling limit of / = 4 works. Besides the three 



The LSND experiment will be ignored throughout this paper. 
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active neutrinos v e , v u , v T on the brane, we must introduce a fourth one, v Sl which is sterile. 
Bearing in mind that oscillation into sterile neutrinos is disfavored, we concentrate on a 
scenario where no induced mixing exists between the sterile and the active brane neutrinos. 
In that case it turns out that the leakage into the bulk comes solely from u s , which is nice 
for two reasons. First, since v s is sterile, there seems to be no reason for it to be confined 
to the brane, so one might expect it to mix more easily with the bulk. Second and more 
importantly, without leakage into the bulk or mixing with u s , the active neutrinos behave as 
if the fifth dimension were absent. That increases the chance for the model to work because 
we know that 4D analyses are phenomenologically successful. Nevertheless, we must realize 
that the presence of the fifth dimension is indispensable for generating the induced mixing 
between the active neutrinos. 

We will show that this decoupling scenario can be obtained by choosing m 4 , the largest of 
the Majorana masses, to be much larger than the other 2/ — 2 = 6 dimensionless parameters. 
Six is just the right number of parameters to fit the three mixing angles and the three 
dimensionless neutrino masses M a R, whatever the radius R is. Nevertheless, there is no 
guarantee that we can get a fit, because the resulting parameters must be ordered in a 
special way to be discussed later. It is therefore encouraging that a fit to experimental data 
is actually possible, with the requisite ordering of the parameters maintained. 

As a consequence of this ordering, we will show that a very interesting connection can be 
established between the smallness of the reactor angle #13 and the smallness of the mass-gap 
ratio AMg 0lar / AM^ tmospheric . 

After reviewing the known result of the strong coupling limit || in Sec. 2 for any /, a new 
parameterization is introduced in Sec. 3, which is needed for the phenomenology application 
in Sec 4. 
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II. THE MODEL AND ITS STRONG COUPLING LIMIT 



The model to be considered here consist of / brane neutrinos, with Majorana masses 
rrii/R, and a single massless bulk neutrino which is coupled to the brane neutrinos through 
Dirac masses di/R, with R being the radius of the extra dimension. In units of 1/R, its 
mass matrix is given by HQ 
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in which the rows and columns are labelled respectively by the left-handed and right-handed 
neutrinos. The first / rows and columns are the brane neutrinos, and the rest are the nth 
mode of the bulk neutrino. The 2f parameters di and m 8 are real and arbitrary, but for 
technical simplicity we will assume the m,;'s to be non-negative, and not equal to an integer 
or a half-integer. 

The eigenvalues A of this matrix satisfies the characteristic equation 

1 



7T 



tan(7rA) = d r(A), 



(2) 



where d 2 = J2i=i d 2 , so e, = di/d obeys the constraint 



/ 

i=i 



1. 



(3) 
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The function r(A) is given by 



r < A > = £rfe < 4) 



In the absence of coupling, the bulk eigenvalues are integers and the brane eigenvalues 
are located at A = m 8 . Accordingly we label the flavor states by an index % or n, where 
i — 1, ■ ■ ■ , / label the brane states and n e Z label the KK tower of bulk states. 

In the strong-coupling limit] 2 ], when d ^> m ; and 1, the characteristic equation (Q) is sat- 
isfied when either tan(7rA) = oo, or r(A) = 0. The former produces half-integer eigenvalues 
A n = n + ~, n G Z. The latter produces / — 1 'isolated' eigenvalues A a (1 < a < / — 1), 
which are the roots of that equation?]. Since r(A) approaches +oo for A = m;+ and — oo for 
A = m»— , it is not difficult to see that the isolated eigenvalues must always occur between 
pairs of m^'s. In other words, 

< mi < Ai < m 2 < A2 < • ■ ■ < A/_i < rrif. (5) 

It is important to realize that in the strong coupling limit, flavor states are so thoroughly 
mixed up that it is impossible to tell whether a mass eigenstate is more brane like or bulk 
like, much less which flavor brane neutrino it should be identified with. In particular, it is 



2 Strong coupling is expected if R~ l is the only additional energy scale that is relevant in the 
problem. This may be the case if there is only one large extra dimension, because then the 
gravitational energy scale is still very high and presumably irrelevant. Neutrino mass is given 
by the formula M = Kvj\[R Q, where v is the Higgs expectation value and k is the brane-bulk 
coupling constant in the Lagrangian. k has a dimension of \^R, and therefore can be expected to 
be of that order if there is no other relevant dimensional quantity around in the bulk. Our coupling 
constant di = MR = kv^/R ~ vR is then much larger than 1, if R is of a sub-milliliter size. 

3 There is a slight change of notation from Ref. [Q]. A2,---,A/ in Ref. [Q] have been renamed 
Ai, • • • , Xf—i- Accordingly the first column of V in eq. (29) of Ref. now appears at the last 
column. 



not true that the states with half-integer eigenvalues are bulk like and the isolated ones are 
brane like. After all, there are only / — 1 isolated eigenvalues and not /. 

The transition amplitude Aij(r) from a brane flavor neutrino j of energy E to a brane 
flavor neutrino i, after it has traversed a distance L = 2EtR 2 , is given by 



Aij{r) 



E u ^ u . 



(6) 



where the infinite- dimensional mixing matrix Ui\ is given by the flavor-i component of the 
normalized mass eigenstate with eigenvalue A. The sum over A n = n + 1 can be computed in 
the strong-coupling limit, resulting in a term proportional to a function g(x), with g(0) = 1 
and g(x) — > (1 — i)/y2nx for large \x\. The final result is 0, 



/-i 



Mr) = E V* a V ja e-^ + V* f V jf g(K 2 r), 



(7) 



a=l 



where K 2 = d 2 (l + d 2 n 2 ) > 1, V ia = Ui\ a , and Vif = e 2 . Explicitly, the / x / induced 
mixing matrix is given by 

e2aW v / si e 2 x 2 2/v/s2 ••• e 2 x 2 j-i/ \/sJ~[ e 2 



V 



\e f x f x/^s{ e f x f2 /^ 
where X{ a = 1/(A Q — m;) and 



(8) 



E 



(9) 



^ (A Q -mi) 2 ' 

The matrix 1/ can be shown to be unitary; actually real orthogonal because and m; are 
assumed to be real. If we replace the function g(K 2 r) in (|7|) by exp(— i\ 2 r), for some real 



number A/, then (|7|) is the usual 4D transition-amplitude formula from flavor neutrino 



n ; 



to rii, through the mass eigenstates n&, with a mixing matrixR V: rij = J2k Vjk"n>k- Moreover, 



4 This relation is usually written as V{ = Ylj=i UijVj (i = e, fi, r) for three flavors. We use this 
unconventional notation to avoid later confusion, because rii turns out to be a permutation of fj, 
and hj a permutation of hj, so U is obtained from V by permutating rows and columns. 



if we choose A/ ^> X a , for 1 < a < / — 1, then for any detector with finite resolution, the 
last term exp(— iX'jr) averages to zero except for tiny r's, just like g(K 2 r). In other words, 
this 5D model in the strong coupling limit resembles a 4D theory with a mixing matrix V 
and a mass pattern (/ — 1) + 1, with a huge mass gap between the first / — 1 neutrinos 
n a and the last one, hf. Since hf represents the total effect of the tower of half-integral 
eigenvalues, rather than oscillating through it, it will be shown immediately below that a 
brane neutrino disappears into it. In other words, we should probably think of A^ as having 
a negative imaginary part. 

Unlike the 4D case, the total probability for an initial flux of flavor j to remain on the 
brane is not 1. There is a leakage into the bulk, equal to 

1 " E \AAr)? = e]{l - \g(K 2 r)\ 2 } ~ e). (10) 
i=i 

The leakage starts out to be at r = 0, but in the strong-coupling limit it reaches its 
asymptotic value e 2 almost instantaneously. This leakage is given by the square of the 
matrix elements in the last column of V. Since J2j = 1 , the total leakage into the bulk 
from all brane neutrinos is 1, independent of /. 

III. A SIMPLE PARAMETERIZATION 

It is fairly complicated to compute the elements of V in (|8|). Given the parameters and 
rrii, we must first find the eigenvalues X a by solving r(A a ) = 0. Then we have to compute 
Xi a and s a , and finally the matrix elements of V. In this section, we discuss a much simpler 
parameterization of V which enables its elements to be computed without having to solve 
any equation. This new parameterization is needed for the phenomenological discussion in 
the next section. 

Instead of the 2f — 1 parameters and m^, we use as independent parameters r/ij and 
X a . Since X a are the zeros of the function r(A) in (^), they must obey the ordering in (^j), 
but otherwise these parameters are free to vary. To show that we must be able to compute 
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e\ from these new parameters. To that end note that the function r(A) is a meromorphic 
function with zeros at X a and simple poles at rrii, and it goes like 1/A for large |A|. Hence 
it is equal to 

n/j(A - x fi ) 



r(A) 



(11) 



Unless otherwise specified, all products over j are taken from 1 to /, and all products over 
(3 are taken from 1 to / — 1. The quantity ef, being the residue at the simple pole m$, is 
then equal to 

2 _ Ufsim ~ X/s) 



.,(///, - ///,)' 

With the ordering (|5]), ef obtained from ([X2|) are automatically non-negative. 
We can now calculate from (|) to get 

_/ -.sg Ylp^aiXa - Xp) 

UjiK-rrij) 



(12) 



(13) 



which is always positive if (§) is obeyed. Since Xi a = l/(X a — rrii), the matrix elements of V 
become 



1 



A a - rrii \ 



Yipijni - Ag) IIj(Ac 



17ln 



2 _ II - V) 



II ; , ,("'/ - m iy 



(14) 



The index % runs from 1 to /, and the index a runs from 1 to / — 1. It is easy to see that 
the sign of V^- is positive unless j < i, when it is negative. With the sign taken care of, the 
magnitude of the matrix elements can be written in the more symmetric form 



\Via\ 

Vif 



\ 



rrin 



2 _ Upim ~ Xp) 



(15) 



Note that V is scale invariant, in that V is unchanged if we multiply every parameter rrii, X 
by the same scale factor. 
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The matrix elements of V in (|15D appear deceptively symmetrical, but they are actually 
not because of (|5|). For example, there is no way to let m 3 — m 2 — > without also letting 
A2 — rri2 and 777.3 — A2 ^ 0. Similarly, there is no way to make 777,3 large without also making 
A3 and 7774 large, but we can make 7774 large without affecting other parameters. 

In view of (|5|), all quantities can be expressed in terms of the 2(/ — 1) consecutive 
distances 



Pa = A a - 777 a , (1 < a < f - 1) 
q a = m a+ i — X a . 



(16) 



Since V^- is scale invariant, there are actually only 2f — 3 independent parameters needed 
to specify V. 

From (|T5|) one sees that every matrix element of V has an equal number of m/s in 
the numerator as the denominator, except V a f and Vf a , which have one more 777/ in the 
denominator than the numerator. Hence if rrif 3> m a , X a , the induced mixing matrix V will 
take on the form 

/ * * • ■ ■ * \ 

* * • • • * 

, (17) 

* * • ■ ■ * 

V •■■ 1 j 

where the asterisks indicate elements that are generally not zero. In that case, there is no 
induced mixing between the first / — 1 flavor neutrinos and the last. There is also no leakage 
of n a into the bulk, for 1 < a < f — 1 . 



IV. PHENOMENOLOGY 

The neutrino mass corresponding to eigenvalue A a is M a = X a /R. In order to have two 
different AM 2 's to explain both atmospheric and solar neutrino oscillations, we need three 
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distinct A Q 's, hence / = 4. Thus in additional to the three brane neutrinos v e , u^, v T which 
are active, we are forced to have one more brane neutrino v s that is sterile. 

In the usual 4D analysis, a sterile neutrino v s is introduced to get us three independent 
AM 2, s, needed to explain solar, atmospheric, as well as LSND oscillations. In the present 
context in five dimensions, leakage into the bulk is simulated by an infinite AM 2 , so we need 
four brane neutrinos just to have two finite AM 2 's. 

Experimentally, oscillation from solar and atmospheric neutrinos into the sterile ones are 
disfavored. Since v s must be present in the theory, this requirement can be fulfilled if v s 
decouples from the active neutrinos, as in fll7l). All we have to do is to let = u s , and 
choose m 4 to be far greater than X a and m a . 

Since V is scale invariant, it is more useful to express it in scale-invariant parameters 



f = Pi/(?i +P2 + q2 + Pz) 

« = +P2 + 92 +P3) 

P = P2I (9l + P2 + 92 + p 3 ) 
7 = 92/ (91 + P2 + 92 + ps) 

5 = P3/(9l +P2 + Q2+P3)- 



(18) 



Only four of these are independent because 



01 + + 7 + 5=1. 



(19) 



See Fig. 1 for a summary of the parameters used. 
The exact expressions for the asterisks in fll?D are: 



(a + £ + /?)(! + Q(7 + a + /?)a 
\ (a + £)(7 + a + £ + /3)(/3 + a) 
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01 + 07/3 



\ (a + 0(7 + a + £ + /?)(5 + 7)(/3 + a) 



13 



Z(a + S + 0)5(i + + 5) 
\ ( a + 0(7 + a + £ + /W + 7) 
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/% + /? + *)(7 + <* + /?)£ 
g(7 + g + dyyja + £ + P) 

(a + 0(7 + /?)(5 + 7) 




7<5a£ 

( 7 + a + e + /5)(7 + /3)(/3 + a) 



(7 + a + + e + (5) 

(7 + « + ^ + /3)( 7 + /3)(5 + 7)(/5 + a) 

( 7 + a + /3) 7 (7 + /? + <*)(l + fl 
(7 + a + e + /^)(7 + /5)(^ + 7) 

(20) 

We shall henceforth denote this 3x3 sub-matrix of V as V. 

With phases neglected, there are six measurable quantities, the three neutrino masses 
and the three mixing angles. There are seven parameters in the model, so one of these 
parameters can never be determined. We shall take it to be the radius R of the fifth 
dimension. The other six parameters m a , X a can in principle be determined from the six 
experimental quantities, but there is no guarantee that the fitted parameters will obey the 
ordering relation (||). It is therefore very encouraging that the resulting fit from the known 
experimental data does obey (|). 

As seen in (BO), four parameters are needed to specify the induced mixing matrix V, 
not three. Hence the values of the neutrino masses M a may affect the values of the mixing 
angles. This turns out to be the case in a particularly relevant situation which we shall 
discuss later. 

The fitting of experimental data is not as straight-forward as it might be, because 
we do not know how V is linked to the experimental mixing matrix U, defined by 
2/j = J2j=i UijVj (i = e, fi, r), and parameterized as 



21 



A 



22 
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V, 
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\ 



12 



C12C13 



S12C13 



Sl3 \ 



u 



-S12C23 - C12S23S13 



C12C23 - S12S23S13 



S23C13 



(21) 



V S12S23 — C12C23S13 



C12S23 — S12C23S13 c 23 c 13 / 



As usual, Sij = sin 8ij and Cjj = cos 6^-, and phases are ignored. We shall follow the usual 
convention to take all mixing angles to be between and ir/2, so that both and cy are 
non-negative. 

We do know that the three n's are flavor neutrinos so they must be the same as the 
three i/'s, and that the three n's are mass eigenstates so they must be identical to the three 
i/'s. But, we do not know which is which, is arranged according to increasing Majorana 
masses mj of the uncoupled flavor neutrinos, and hj is arranged according to increasing 
eigenvalues. On the other hand, v\ is fixed by their weak interactions, and Vj is tied to V{ 
through Uij and the experimental mixing angles. To find the correct assignments of rij to z/j 
and hj to Uj, we need to determine what permutations of rows and columns of V are needed 
to bring it into the form of U. As noted earlier, although there are six parameters in V, the 
ordering requirement ([5]) does make the columns and rows of V inequivalent. 

The present experimental situation is as follows. 

Atmospheric mixing is maximal. It gives tan 2 823 = 1, with a square mass difference 



most likely ones: LMA, LOW, and VAC. The exact parameters and the goodness of the 
fit of each depend on the specific analysis, but generally speaking, the LMA solution is the 
most favorable. In global analyses where both rate and spectrum are taken into account, 
the SMA solution is poor so we will not consider it. The parameters do not vary that much 
from analysis to analysis, so for the following discussions, we will adopt the parameters taken 



from Table 2 of Ref. 0. The (tan 2 # 12 , AM s 2 ol ) values for these solutions are: LMA=(0.36, 
5 x 1(T 5 ), VAC=(0.363, 1.4 x 10" 10 ), LOW=(0.69, 1.1 x 1(T 7 ). The unit for AM 2 is (eV) 2 . 



AM 2 
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As for the reactor angle, only an upper bound of tan 2 #13 < 0.04 is knowr|] from the 



CHOOZ experiment [ [TT|j . 

We do not know the values of M a , except for an upper bound of about 2 eV |12[] from 
end point measurements of tritium /?-decay. We will therefore treat Mi as a parameter, 
to be varied from to 2 eV. Using (|^), the values of M 2 and M 3 are then determined 
from the magnitude of the solar and atmospheric gaps. If the solar gap is on top, then 
M 2 = J Ml + (AM 2 ) atm and M 3 = J Ml + (AM 2 ) sol . We will designate this case as G = 1. 



If the solar gap is at the bottom, then M 2 = yjM% + (AM 2 ) so , and M 3 = ^/Mf + (AM 2 ) aim . 
We will designate this case as G = 2. 
Let 



p=(M 2 - Mi)/(M 3 -M 1 ) = a + j3 = l-'y-6 



(22) 



be the ratio of mass differences. Its value depends on Mi, on G, and on the solar- neutrino 
solution. Listed below are some values of p for the solar solutions [LMA, LOW, VAC] at 
various Mi and G: 

[.98995, .99998, 1.00000] (M l = OeV, G = 1) 
[.98001, .99996, 1.00000] (Mi = leV, G = 1) 
[.98000, .99996, 1.00000] (M x = 2eV, G = 1) 
[.14142, .00663, 0.00024] (Mi = OeV, G = 2) 
[.02001, .00004, 0.00000] (Mi = leV, G = 2) 

[.02000,.00005,0.00000] (Mi = 2eV, G = 2). (23) 



We see that p is close to 1 for G = 1, and close to for G = 2. 



5 CHOOZ quotes an upper bound of sin 2 (2#) = 0.1 for large AM 2 , corresponding to a tan 2 9 = 
0.026. In the vicinity of AM 2 = 2.5 x 10~ 3 eV 2 , the best-fitted atmospher ic mass difference, the 
bound is tan 2 9 ~ 0.038. We will be cautious and take tan 2 #i 3 < 0.04. 
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An analytical solution exists if we approximate p to be for G = 2, and 1 for G = 1. This 
is discussed below, separately for G = 2 and G = 1, together with the numerical solutions 
when this approximation is not made. 



A. G = 2 



The mass of rii increases with i. When the solar gap is at the bottom, either (ni, n<i, ^3) = 
(^1,^2,^3), or (ni,n 2 ,n 3 ) = (i> 2 , z>i, z> 3 ). 

When p = 0, it follows from (|22|) that a = /3 = 0, because both a and /3 must be 
non-negative. Using (|19l), we also get 5 = 1 — 7. Defining a/ (3 = a, the matrix V in this 
limit simplifies to 

I ^(1 + 07^/(7 + 0(1 + 0) ^(1 + 07/(7 + 0(1 + a) V / ^ 1 -^/^ + ^\ 





a/(l + a) 



V-,/(l- 7 K/(7 + 0(l + a) -t/(l- 7)^/(7 + 0(1 + a) J 7 (l + 0/(7 + 0/ 



(24) 



The outstanding feature of this matrix is V^23 = 0. Since n 3 = P3 and the atmospheric 
mixing is maximal, the only way that this could happen is for n<i = v e , and V23 = U e s = 

S13 = 0. 

Thus the smallness of S13 is naturally related to the smallness of AM s 2 oj /AM^ m , at least 
in the sense that we are forced to have S13 = if the mass ratio is zero. This remarkable 
feature of the model is a consequence of the ordering (|), which demands both a and (3 to 
be non-negative. It is not the result of an arranged texture. 

To implement the atmospheric mixing condition tan 2 #23 = 1, we need £ = 7/(1 — 27), 
then 



V 



f y/a/2(l + a) ^1/2(1 + a) .Jl/2\ 





1/(1 + a) 



a/(l + a) 



(25) 



\-y/a/2{l + a) -^1/2(1 + 0) v/1/2/ 
To bring into the form of £7, we can multiply the first column of V by —1 (this is just a 
change of sign convention), and change its rows (1,2,3) to (2,1,3). The resulting matrix 
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/ ^1/(1 + a) y/a/(l + a) \ 



U= -y/a/2(l + a) yjl/2(l + a) y/T/2 



(26) 



V y/a/2(l + a) -^1/2(1 + a) y/T/2 J 



should be compared with U in the approximation si 3 = 0, which is 



/ C12 



\ 
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-S12C23 



C12C23 



(27) 



V S12S23 -C12S23 c 23 / 



It shows S12 = Ja/(1 + a) and C12 = J 1/(1 + a). Since a is arbitrary, all solar solutions can 
be accommodated. 

If we interchange the first two columns and also the last two rows of U, we simply 
interchange a with a~ l . Since a is arbitrary, this form is equally allowed. Hence there 
is no way to distinguish the case (n 2 ,n 3 ) = (i/^,i/ t ) and (ni,n 2 ) = (^1,^2) from the case 
(n 2 ,n 3 ) = (wr,*v) and (ni,n 2 ) = (£ 2 ,r>i). 

The actual value of p is given in (|23| ) . It is zero for the VAC solution when Mi is between 
1 and 2 eV, to within the number of digits shown, so the analytical solution discussed above 
is perfectly adequate. No more will be said about it. 

For the other cases, p > 0, numerical solutions are provided in Table I. This Table is 
constructed in the following way. For a given M\ and £, and a given solar solution, the 
three parameters /3, 7, 8 (a = 1 — (3 — 7 — 5) are determined from fitting the three inputs, 
p, tan 2 #23 = 1, and tan 2 612. Then tan 2 #13 is calculated. To test the goodness of the 
approximation p = 0, we also list in the last two columns a = a//3 and 7/(1 — 27) calculated 
from the fitted parameters. When p = 0, the former is equal to tan 2 # 12 , and the latter is 
equal to £. Other than the LMA solution at M\ = eV, where p = .14 is fairly large, the 
agreement is quite good. This shows that the approximation ]rho = is a very reasonable 
one. 

Since tan 2 #13 = when p = 0, one might expect tan 2 6*13 (and therefore also sin 2 #13) to 
be small when p is small. A glance at Table I shows that this is not necessarily the case, 
as we can adjust the parameter £ to yield a tan 2 #13 as large as the CHOOZ bound 0.04, or 



16 



even bigger. To understand what this means let us examine the formula 

- 2 °* = Vis = 1 t 5{l + Jl ^ ( 2 8) 

obtained from (p0|). For small p, both a and (3 are small. As pointed out above, the 
approximations £ ~ 7/(1 — 2j) and tan 2 612 — a = a/ (3 are good. The parameter a is fixed 
by the solar solution, so sin 2 # 13 is of order /3 2 , unless the denominator in (f28[) is also small. 
Now 5 + 7= 1 — a — (3 ~ 1 for small (3 and a, and £ is proportional to 7 for small 7, hence 
the denominator is small and of order (3 2 if and only if 7 is itself of order (3. In that case 
£ ~ 7, and 

sin 2 g 13 ~ ° ( 7 = c /5«l) (29) 

(a + c)(l + c) 

may indeed be large. An examination of Table I shows this is indeed what happens: tan 2 #13 
is relatively large only when 7 is relatively small. 

7 can vary from Otol — a — /3 ~ 1 for small p. Except for very small 7, in most of 
this range it yields a small tan 2 # 13 . So in a statistical sense it is indeed true that a small p 
tends to yield a small tan 2 # 13 , meaning that the smallness of the reactor angle is naturally 
related to the smallness of the solar to atmospheric mass gap ratio. 



B. G = 1 

The main conclusions of the last subsection are essentially unchanged. Here are the 
details. 

When the solar gap is on top, either (n 1; n 2 , n 3 ) = (z> 3 , P2, V\) or (n l5 n 2 , n 3 ) = (i> 3 , v\, v 2 ) 
must be true. 

When p — 1, it follows from ( p2|) that 7 = 5 = because they must both be positive. In 
that case it follows from flT9|) that a = 1 — /3. Defining 5/7 = cf, the matrix 1/ in this limit 
simplifies to be 
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Similar to the case G = 2, since n\ = u s , we should identify V31 = U e3 = si 3 = 0, and hence 
n 3 = v e . The atmospheric condition tan 2 #13 = 1 then requires £ = (1 + d — [3 — d(3)J (—1 — 
d + 2(5 + df3), and turns V into 

/ y/l/(2 + d) y/l/(2 + d) y/d/(2 + d) \ 



V 



(l + d)/(2 + d) x /l/(l + d)(2 + d) x /d/(l + d)(2 + d) 




(31) 



V -yfdToTd) + d) j 

Multiplying the second row and the second column each by —1, and interchanging the first 
and third rows and well as the first and third columns, bring V to the form 
/ y/W + d) y/d/{l + d) \ 

1= - Jd/(l + d)(2 + d) Jl/(l + d)(2 + d) J(l + d)/(2 + d) , (32) 



- A /d/(l + d)(2 + d) y/l/(l + d)(2 + d) y/(l + d)/(2 + d) 

V y/d/(2 + d) -VV(2 + d) y/l/(2 + d) ) 

which we can identify with U of (|27|) to fix the parameter d from the solar angle. 

we can carry out a numerical fit like Table I and obtain similar results. Since in this case 
p ~ 1 is a fairly good approximation even for the Mi = 2 situation of LMA, we will not 
show these fits here. 



V. CONCLUSION 

In theories with extra dimensions where ordinary particles are confined to the three- 
brane, only gravity and sterile neutrino are allowed to roam in the bulk. Moreover, if there 
is only a single large extra dimension, the energy scale of gravity would still be too high to 
produce many gravitons, then it leaves neutrino to be the only probe of extra dimensions 
in high-energy experiments. It is therefore important to investigate whether there are signs 
pointing to a higher dimension in the neutrino data. 
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A large extra dimension may explain the small neutrino mass. Strong coupling between 
brane and bulk neutrinos may also explain large mixing not seen in the quark sector. Un- 
fortunately, it is not easy to realize these ideas, because oscillation into sterile neutrinos 
are greatly disfavored in the solar and atmospheric neutrino data, and the presence of the 
Kaluza-Klein tower of sterile neutrinos from the bulk tends to produce a very complicated 
oscillation pattern not hitherto observed. In this paper we consider a simple and economic 
model where both of these problems are solved. The model consists of a sterile neutrino v s 
in additional to three active neutrinos is e , z/^, v T on the brane, coupled strongly to a single 
common neutrino in the bulk, but not directly among themselves. In the strong coupling 
limit, the complicated oscillation patterns are washed out, the induced couplings between 
the brane neutrinos become large, and the sterile neutrinos on the brane and in the bulk can 
be made to decouple from the active neutrinos. This results in a realistic theory of neutrino 
oscillations, capable of accommodating all existing neutrino data. It also predicts a natural 
connection between the smallness of the reactor angle 9 13 , and the smallness of the solar to 
atmospheric mass gap ratio. 

This research is supported by the Natural Science and Engineering Research Council of 
Canada and FCAR of Quebec. 
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TABLES 

TABLE I. Each of these tables shows the fitted parameters (a, 0, 7, S) for a given M\ and £, 
when G = 2. The resulting value of tan 2 #13 is given. The last two columns should be equal to 
tan 2 O12 and £ respectively in the p = limit. To save space, 1.8 x 10~ 5 , for example, is written as 
1.8(-5). 



LMA. Mi = 2 eV. p = .02. tan 2 12 = 0.36. 



s 


tan 2 #13 


a 





T 


s 


a/0 


7/fl - 27) 


.500 


.0007 


.0054 


.0146 


.2465 


.7335 


.3668 


.4862 


.300 


.0013 


.0054 


.0146 


.1868 


.7932 


.3662 


.2983 


.050 


.0211 


.0052 


.0148 


.0520 


.9280 


.3488 


.0581 


.032 


.0398 


.0050 


.0150 


.0375 


.9425 


.3329 


.0406 


LMA. Mi 


= eV. p = 


.14. tan 2 # 12 


= 0.36. 










S 


tan 2 #13 


a 





7 


s 


a/0 


7/(1 - 27) 


100 


.0043 


.0415 


.0999 


.4106 


.4480 


.4150 


2.297 


10 


.0050 


.0414 


.1000 


.3937 


.4648 


.4142 


1.853 


1 


.0139 


.0408 


.1006 


.2846 


.5740 


.4054 


0.661 


0.34 


.0391 


.0390 


.1024 


.1901 


.6685 


.3807 


0.307 


LOW. Mi 


= 2 eV. p = 


.00005. tan 2 


Q X2 = 0.69. 












tan 2 #13 


a 





7 


5 


a/0 


7/(1 - 2 7 ) 


.5000 


4.4(-9) 


1.8(-5) 


2.7(-5) 


.2500 


.7500 


.6900 


.4999 


.0100 


4.9(-6) 


1.8(-5) 


2.7(-5) 


.0098 


.9901 


.6900 


.0100 


.0010 


.0005 


1.8(-5) 


2.7(-5) 


.0010 


.9989 


.6898 


.0010 


.0001 


.0313 


1.8(-5) 


2.7(-5) 


.0001 


.9998 


.6738 


.0001 
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LOW. M 1 = eV. p = .00663. tan 2 # 12 = 0.69. 



e 


tan 2 #i3 


a 





7 


5 


a/0 


7/(1 " 2 7 ) 


i 


4.2(-5) 


.0027 


.0039 


.3302 


.6632 


.6946 


.9723 


.100 


.0012 


.0027 


.0039 


.0833 


.9101 


.6940 


.0999 


.030 


.0099 


.0027 


.0039 


.0291 


.9643 


.6894 


.0309 


.013 


.0398 


.0027 


.0040 


.0137 


.9797 


.6739 


.0141 


VAC. 


Mi = eV. p = 


.00024. tan 2 


12 = 0.363. 












tan 2 6*13 


a 





7 


5 


a/0 


7/(1 " 2 7 ) 


.1000 


1.3(-6) 


6.3(-5) 


.0002 


.0834 


.9164 


.3631 


.1001 


.0100 


.0001 


6.3(-5) 


.0002 


.0099 


.9899 


.3630 


.0101 


.0010 


.0080 


6.2(-5) 


.0002 


.0011 


.9987 


.3562 


.0011 


.0004 


.0400 


5.9(-5) 


.0002 


.0005 


.9993 


.3292 


.0005 
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FIGURES 
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FIG. 1. A summary of the parameters in the strong-coupling limit for / = 4. rrn are the 
Majorana masses and X a are the isolated eigenvalues. These parameters must be ordered as 
indicated. Shown in the figures are also the distances p a ,q a as well as the ratios £, a,/3, 7, and 5. 
The common distance with respect to which the ratios are taken is A3 — Ai. 
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